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Abstract 

Inclusive J/ip production in photon-photon collisions has been observed at 
LEP II beam energies. A clear signal from the reaction 77 — > J/ip + X 
is seen. The number of observed NlJ/ip — > yU^yU") events is 36 ± 7 for an 
integrated luminosity of 617 pb7^ yielding a cross-section of uiJ/ip + X) = 
45 ± 9(stat) ± 17(syst) pb. Based on a study of the event shapes of different 
types of 77 processes in the PYTHIA program, we conclude that (74 ± 22)% 
of the observed J/ip events are due to 'resolved' photons, the dominant contri- 
bution of which is most probably due to the gluon content of the photon. 
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1 Introduction 

An important component of the e'^e' collisions at LEP II energy is the two-photon 
fusion process. It has been pointed out that two-photon production of inclusive J/ip 
mesons: 

e"*" + e+ + e~ + 7i + 72, (1) 

7,+7, ^ J/V^ + X. (2) 

is a sensitive channel for investigating the gluon distribution in the photon [1]. 

There are two important processes leading to inclusive J/ip production. The corre- 
sponding typical diagrams are given in Fig. la-b. Less important diagrams are not 
considered here. The first process is described by the vector-meson dominance (VDM) 
model [2]: 

7i -> c + c, 72 9 + ^ and 1, + % ^ J/^ + Q + (3) 

as shown in Fig. la. The vertices for 7^ and 72 are connected by Pomeron exchange or 
diffractive dissociation of photons. The final-state parton pairs c + c and q + q are both 
in the state of J^'" = 1 , which means that the latter is dominated by the low-mass 
vector mesons p°, u and but a more general inclusive hadronization of the partons may 
also be important. 

The second process is described, for example, by the colour-octet model [3]. It proceeds 
through the so-called 'resolved' contribution of the photons, in which the intermediate 
photons are 'resolved' into their constituent partons: 

7i+ 9^ ^ c + c, ^2+ 9^ ^ Q + Q, and 7i + ^ J/i^ + Q + ^, (4) 
as shown in Fig. lb. 



a) b) 




Figure 1: Inclusive J/tp production in 77 processes: a) through vector-meson dominance, 
b) via gluon content of the photon, i.e. 'resolved' contributions. 

It is seen that this process requires production of a 'resolved' gluon ((7^) from both 
photons. Thus, this production mechanism provides a sensitive probe of the gluon content 
of the photon. 

The purpose of this letter is to report the observation of inclusive J/ip production 
from the two-photon fusion process, to give its production characteristics along with the 
cross-section and finally to assess the relative importance of the production processes 
discussed above. Section 2 describes the selection criteria for the event sample collected 
for this study. The measurement of inclusive J/ip production in the fi~^fi' channel and 
its interpretation in terms of diffractive and resolved processes is presented in section 3 
followed by a summary and conclusions. 
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2 Experimental Procedure 

The analysis presented here is based on the data taken with the DELPHI detector [4,5] 
during the years 1996-2000, excluding the part of the data collected in the last period of 
2000, when one of the Time Projection Chamber (TPC) sectors was not in operation. The 
centre-of-mass energies \fs for LEP ranged from 161 to 207 GeV. The total integrated 
luminosity used in the analysis is 617 pb~^. 

The charged particle tracks were measured in the 1.2 T magnetic field by a set of 
tracking detectors including the micro Vertex Detector (VD), the Inner Detector (ID), 
the TPC, the Outer Detector (OD) and the Forward/Backward Chambers FCA and 
FCB. The following selection criteria were applied: 
(a) particle momentum p > 200 Me\/c; 
(6) relative momentum error of a track Ap/p < 100%; 

(c) impact parameter of a track, transverse to the beam axis < 3 cm; 

(d) impact parameter of a track, along the beam axis < 7 cm; 

(e) polar angle of a track, with respect to the beam axis 10° < 9 < 170°; 
(/) track length > 30 cm. 

The neutral particles (7, 7r°, K^, n) were selected by demanding that the calorimetric 
information, not associated with charged particle tracks, satisfies the following cuts: 

[g) i?(neutral) > 0.2 GeV for the electromagnetic showers, unambiguously identified 
as photons; 

{h) i?(neutral) > 0.5 GeV for all the other showers; 

(i) polar angle of neutral particle tracks, with respect to the beam axis 
10° < ^ < 170°. 

In order to ensure a very high trigger efficiency, the selected events were required to 
satisfy at least one of the following sets of criteria: 

(ji) one or more charged particle tracks in the barrel region (40° < 6 < 140°) with 

Pt > 1.2 Gey/c, is found; 
(j2) one or more neutral particle tracks in the Forward ElectroMagnetic Calorimeter 
(FEMC) ( 10° < ^ < 36° and 144° <$ < 170°) with energy greater than 10 GeV, 
is found; 

(js) the total sum of charged particle tracks in the barrel with pt > 1 GcV/c, of charged 
particle tracks in the forward region (10° < ^ < 40° or 140° < 9 < 170°) with 
Pt > 2 GeV/c and of neutral particle tracks in the FEMC with E > 7 GeV, is 
greater than one; 

(j^) the total sum of charged particle tracks in the barrel with pt > 0.5 GeV/c, of 
charged particle tracks in the forward region with pt > 1 GeV/c and of neutral 
particle tracks in the FEMC with E > 5 GeV, is greater than four. 

The trigger efficiency for the events which passed the above requirements is bigger than 

98%. 

The hadronic two-photon events are characterized by a low visible invariant mass. 
Consequently, the following additional cuts were applied: 

(k) the visible invariant mass, Wyis, calculated from the four- momentum vectors of 

the measured charged and neutral particle tracks, is less than 35 GeV/c^; 
(/) the number of charged particle tracks N^h satisfies 4 < A^ch ^ 30; 
(m) the sum of the transverse energy components with respect to the beam direction 
for all charged particle tracks (Yl \JVt ^ ^l) is greater than 3 GeV. 
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The comparison of the W^is distributions, after the cuts on N^h and ^ -E'™, both for the 
data and the events simulated by PYTHIA, shows (Fig. 2) that the cut W^is < 35 GeV/c^ 
rejects the major part of the non-two-photon events. 




Figure 2: W^is distributions for the LEP II DELPHI data, for the simulated 
77 — hadrons, e~^e~ — > Z^'-f, e~^e~ — > W^W~ and the sum of all above Monte-Carlo 
contributions 

A total of Nt = 274 510 events remain in the data sample after applying all these cuts. 
The main background comes from the process e~^e~ Z'y and amounts to ~1.2% of the 
selected 77 events. The background from the e'^e" W^W~ is negligible, as seen in 
Fig. 2. 

J/ip candidates have been selected using the fi^fi^ decay channel. For the muon pair 
selection, the following criteria were imposed: 

(n) at least two charged particle tracks, with zero net charge, should be accepted by 
the standard DELPHI muon-tagging algorithm [5], or be identified as muons by 
the hadronic calorimeter; 
(o) the tracks should not come from any reconstructed secondary vertex or be 
identified as a kaon, proton or electron by the standard DELPHI identification 
packages. 

3 Inclusive J/if) Production 

In this section, we first determine the inclusive J/i/j production in the fi~^fi~ chan- 
nel. Then we interpret it in terms of diffractive and resolved processes by fitting the 
experimental p'^{J/ip) distribution to the PYTHIA predictions for both processes. This 
allows to deduce the cross-section for inclusive J/ip production, taking into account the 
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7 7 + ^ the J/il] — >• efficiencies. As the ffist set of efficiencies is model- 

dependent, we also give the 'visible' cross-section in which only the detector efficiency for 
J/%1) — > decay is considered. We finally present the production characteristics 

together with the PYTHIA predictions. 

In Fig. 3 we give the invariant mass distribution for identified fi'^^~ pairs, selected 
as outlined in the previous section. The J/ip signal shows up over little background. A 
least squares fit to the M{iJ,'^iJ,~) distribution with a Gaussian for the signal and a second 
order polynomial for the background gives the following results: 



J/ip mass: M 
J/ip width: (7(obs) 



3119±8 MeY/c^, 
35±7 MeV/c^. 




3.2 3,4 3,6 3,8 

M(/iV). GeV/c 



Figure 3: M{fi^fi ) distribution from the LEP 11 DELPHI data. The solid curve corre- 
sponds to a Gaussian fit over a second order polynomial background. 



The observed width of the peak is consistent within errors with the invariant mass 
resolution of a pair of charged particle tracks in the mass region around 3Ge\/c^. The 
number of observed events from the fit is: 

N{J/iIj) = 36±7 events, 

over a background of about 11 events. 

If we take the L3 result [6] for the beauty cross-section from 77 events and the PDG 
value [7] for the branching ratio of beauty hadrons to J/ip, the expected number of 
J/ip —>■ ^x^^JL~ from beauty hadrons is 2.1 ± 0.6. The backgrounds from the processes 
e+ -I- e~ — > Z -I- 7 — >■ J/ip + X and 7 + 7 — Xc2 — ^ J/i' -|- tt"*" -|- tt" + 7r° are less than 0.20 and 
0.30 event respectively. According to the selection criteria the system X contains at least 
two charged particle tracks, hence wc do not consider such sources of J/ip production as 
7 4- 7 — Xc2 — ^ + 7- We checked that in the four-prong events with J/%1) A*"*"/^" 
candidates there are no photon conversions. 
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We used the PYTHIA 6.156 generator [8] to estimate the efficiency. The generated 
events were passed through the simulation package of the DELPHI detector [5] and then 
processed with the same reconstruction and analysis programs as the real data. There is 
a substantial fraction of PYTHIA events where J/tl) mesons are produced just as a simple 
fusion of two photons because there is not enough phase space to produce additional 
particles. We do not use such events. The process where both photons are VDM photons 
we will call 'diffractive' and the process without VDM photons we will call 'resolved'. 



> 



Q. 

< 



DELPHI 

Diffractive 
Resolved 




Figure 4: p^iJ/ip) distribution from the LEP II DELPHI data, shown as points with 
error bars. The histogram is a combination of the normalized 'resolved' and 'diffractive' 
processes from PYTHIA (see text). 



A set of the J/t/j production characteristics is exhibited in Figs. 4 to 8. For each bin of 
every distribution shown, wc have examined the M{/j^fjr) spectrum and then fitted with 
a Gaussian and a second order polynomial background, to get the number of signal events 
per bin. This number is then renormalized for each distribution, so that the total sum is 
always equal to 36 events. Hence, Fig. 4 to Fig. 8 are background subtracted distributions. 

Fig. 4 shows the p^^{J/ip) distribution. As expected, the PYTHIA Monte Carlo pre- 
diction for the p^iJ/ip) distribution is more sharply peaked near zero for the 'diffractive' 
events (see Fig. la) than for the 'resolved' events (see Fig. lb). We fitted the experimental 
PtiJ/i^) distribution as a function of the two categories of PYTHIA events: 



dA^ 

T 



dp: 



/ 



dA^ 



dp: 



+ (!-/)• 



dA^ 



Diffractive 



dp: 



(5) 



Resolved 



which gives / = (26 ±22)% (PYTHIA distributions in Fig. 4 are normalized to the data). 
The PYTHIA study tells us that the experimental efficiencies are very different for the 
two categories: 



e( diffractive) 
e(resolved) 



(0.98 ± 0.04)%, 
(3.87 ±0.09)%. 



(6) 



According to PYTHIA, about one-half of all the 77 events with J/ip — > yU+yU are produced 
with the charged particle tracks at polar angles below 10 degrees, so that they are invisible 
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to the DELPHI detector. The individual efficiencies as a function of are given in Fig. 5. 
Some insight may be gained into these efficiencies if they are broken down into products 
of two factors, as follows: 



e(diffi^active) 
e (resolved) 



( diffract ive) x e 
(resolved) x e 



^77 



_ (diffractive), 
(resolved), 



(7) 



where e-yy is the efficiency for the process 77 —> J/i^ + X and ^^^_^^+^- is that for 
J/%1) — > //"^/x". As expected, the latter is relatively process- independent: 



. (diffractive) = (37.0 ± 1.5)%, 
(resolved) = (32.1 ± 0.7)%. 



It is clear, therefore, that the difference in efficiency in (6) is mostly due to e^-y. 
highly pro cess- dependent and hence model-dependent. 
The overall experimental efficiency is: 



/ 



+ 



1-/ 



e(diffractive) e(resolved) ' 



(8) 
This is 

(9) 



which gives e = (2.19_o;59)%. Under the assumption that PYTHIA captures the 
kinematical features of the resolved and diffractive processes, but not their absolute cross- 
sections, the cross-section for inclusive J/tl) production is: 



a = N{J/ij) -{Br-C- e)'^ = 45 ± 9(stat) ± 17(syst) pb. 



(10) 



where Br = (5.88±0.10)% is the branching ratio for J/-ip ^ ji'^ji' [7] and C = 617pb~^ is 
the total integrated luminosity. The systematic uncertainties include both the efficiency 
(9) and the branching ratio contributions but not those inherent to the PYTHIA program. 



o 
c 



0.06 



0.02 



DELPHI 



Resolved 



A— Diffractive 



23456789 10 
2 l'n^\l / 



P/, (GeV/c)^ 

Figure 5: Efficiencies for resolved and diffractive processes as functions of p^. 

Because of the model-dependent aspect of this analysis (see, for example, the effi- 
ciencies given in (6)), it is of interest to quote the 'visible' cross-section. Substituting 



7 



(diffractive) and ^jj^^^+^j^- (resolved) for e( diffract ive) and e(resolved) respec- 
tively in (9), the 'visible' cross-section can be calculated; it is: 

(7vis = 3.0 ± 0.6(stat) ± O.l(syst) pb. (11) 

The main source of systematic uncertainty comes from the determination of the relative 
fractions of resolved and diffractive events which have different efficiencies (8). Following 
the same argument, we also give the 'visible' production rate (n) for J/%1) production: 



(n) =N{J/ij) ■ [NfBr 



(6.7 ± 1.3(stat) ± 0.3(syst)) x 10' 



(12) 



where A^^^ is the data sample for the 77 selection as given in the previous section. 

The rapidity distribution in the laboratory system for the J/ijj mesons is shown in 
Fig. 6. The PYTHIA events have been combined using the same fraction / found in 
(5) and then normalized to the observed number of events in < \y\ < 2.0. The same 
techniques have been used to compare the experimental distributions of M{J/iIj + X), 
M{X), the charged and total multiplicities [Nch{X) and A^tot(^)], in Figs. 7a-d. There 
is fair agreement within statistics between the shapes of our measured distributions and 
the PYTHIA predictions (using the best fit as found in (5) for the relative content of 
diffractive and resolved events and renormahzing the PYTHIA prediction to the number 
of observed events) . 



DELPHI 



Diffractive 
Resolved 




Figure 6: \y\ distribution for J/tjj mesons from the LEP II DELPHI data, shown as 
points with error bars. The histogram is a combination of the normalized 'resolved' and 
'diffractive' processes from PYTHIA (see text). 

The acceptance-corrected distributions in cos^, where 6 is the helicity angle of fi"^ 
in the rest frame of J/ip — > A*^/^", are shown in Figs.8a-c, along with the results of 
a fit to the form {1 + a cos^O). The fitted parameters are a — —0.9 ± 0.6 for the 
total sample (a), a = -1.8 ± 0.5 for (j/^) < 1.0 (GeV/c) (b) and a = 0.7 ± 1.3 for 
p'^{J/ip) > 1.0 {GeV/cY (c). These results indicate that the J/ip mesons are produced with 
little polarization at high p^iJ/ip), where the main contribution comes from the resolved 
processes. 
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Multiplicity of X 



Figure 7: Visible distributions of M^J/ip + X) (a), M(X) (b), charged (c) and total 
(d) multiplicities of the X system. Each histogram is a combination of the normalized 
'resolved' and 'diffractive' processes from PYTHIA (see text). 




Figure 8: Acceptance-corrected distributions in cos 9 where 6 is the helicity angle of /i^ 
in the rest frame of J/if) — > /^"""/i". The figures (a-c) correspond to the total sample (a), 
the subsamples with (j/^) < i.o (GeV/c)^ (b) and > i.o {GeY/cf (c). 
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4 Conclusions 

We have studied the inclusive J/-^ production from 77 coUisions. The data have been 
taken by the DELPHI Collaboration during the LEP II phase, i.e. ^/s of the LEP machine 
ranged from 161 to 207 GeV. A clear signal from the reaction 77 J/ip + X is seen. 

The inclusive cross-section is estimated to be a{J/ip + X) =45 ± 9(stat) ± 17(syst) pb. 
Based on a study of the distribution of different types of 77 processes in the PYTHIA 
program, we conclude that some (74 ± 22 )% of the observed J/ip events are due to the 
'resolved' photons, the dominant contribution of which should correspond to the gluon 
content of the photon [3] . 

The distributions in pl,{J/ip), y and cos^ (for in the rest frame of J/ijj 
are presented. In addition, a study is given of the characteristics of the system X. 
All distributions appear to be well reproduced within statistics by the normalized 
combination of the fitted 'resolved' and 'diffr active' contributions. 

Acknowledgements 

We are greatly indebted to our technical collaborators, to the members of the CERN- 
SL Division for the excellent performance of the LEP collider, and to the funding agencies 
for their support in building and operating the DELPHI detector. 
We acknowledge in particular the support of 

Austrian Federal Ministry of Education, Science and Culture, GZ 616.364/2-III/2a/98, 
FNRS-FWO, Flanders Institute to encourage scientific and technological research in the 
industry (IWT), Federal Office for Scientific, Technical and Cultural affairs (OSTC), 
Belgium, 

FINEP, CNPq, CAPES, FUJB and FAPERJ, Brazil, 
Czech Ministry of Industry and Trade, GA CR 202/99/1362, 
Commission of the European Communities (DC XII), 
Direction des Sciences de la Matiere, CEA, France, 

Bundesministerium fiir Bildung, Wissenschaft, Forschung und Technologic, Germany, 
General Secretariat for Research and Technology, Greece, 

National Science Foundation (NWO) and Foundation for Research on Matter (FOM), 
The Netherlands, 
Norwegian Research Council, 

State Committee for Scientific Research, Poland, SPUB-M/CERN/PO3/DZ296/2000, 

SPUB-M/CERN/PO3/DZ297/2000 and 2P03B 104 19 and 2P03B 69 23(2002-2004) 

JNICT- Junta Nacional de Investigagao Cienti'fica e Tecnologica, Portugal, 

Vedecka grantova agcntura MS SR, Slovakia, Nr. 95/5195/134, 

Ministry of Science and Technology of the Republic of Slovenia, 

CICYT, Spain, AEN99-0950 and AEN99-0761, 

The Swedish Natural Science Research Council, 

Particle Physics and Astronomy Research Council, UK, 

Department of Energy, USA, DE-FG02-01ER41155, 

EEC RTN contract HPRN-CT-00292-2002. 

We are indebted to Dr. T. Sjostrand for his help with the diagrams included in this paper 
and for his comments on PYTHIA. We thank M. Klasen for his help with the diagrams 
important for inclusive J/ip production. S.-U. Chung is grateful for the warm hospitality 
extended to him by the CERN staff during his sabbatical year in the EP Division. 



10 



References 

[1] 'Physics at LEP2,' edited by G. Altarelli, T. Sjostrand and F. Zwirner, 
CERN96-01 (Vol. 1), Feb 1996. 

Sec Section 7 (heavy-quark physics), p. 330, in the chapter on 77 physics. 
[2] J.J. Sakurai and D. Schildknecht, Phys. Lett. B40 (1979) 121; 

I. F. Ginzburg and V.G. Serbo, Phys. Lett. B109 (1982) 231. 
[3] R. M. Godbole et al, Phys. Rev. D65 (2002) 074003; 

M. Klasen et al, NucL Phys. B609 (2001) 518; 

B. Naroska, NucL Phys. B (Proc. SuppL) 82 (2000) 187; 

H. Jung, G. A. Schuler and J. Terron, Int. J. Mod. Phys. A7 (1992) 7955; 

E. L. Berger and D. Jones, Phys. Rev. D23 (1981) 1521. 
[4] P. Aarnio et al, DELPHI CoUab., NucL Inst. Meth. A303 (1991) 233. 
[5] P. Abreu et al, DELPHI Collab., NucL Inst. Meth. A378 (1996) 57. 
[6] R. Acciarri et al, L3 Collab., Phys. Lett. B503 (2001) 10. 
[7] Particle Data Group, K.Hagiwara et al., Phys. Rev. D66, 010001 (2002). 
[8] T. Sjostrand et al, Comput. Phys. Comm. 135 (2001) 238. 



